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SATELLITE PASSIVE REMOTE SENSING OF OFF-SHORE POLLUTANTS 


INTRODUCTION 

In this report we discuss satellite detection and monitoring of off-shore 
dumped pollutants, other than oil. Following EPA guidelines, we divide 
the pollutants into four categories: 

e acid waste 
9 industrial waste 

• sewage sludge 

• dredge material 

Off-shore dumps consist of material that sink, or at the very least do 
not remain on the water surface. Of the three portions of the spectrum, 
visible, infrared and microwave, only the first that penetrates the water 
can sense waste concentration. The major problem in the visible band is 
to get around weather and atmospheric haze. 

The results of the analyses in this report confirm the intuitive notion 
that all three satellite based sensors have the required sensitivity to 
do the job, but only the visible has sufficient spatial resolution. How- 
ever, none of the sensing techniques allow a clean cut extraction of the 
pollutant signature from the background. We assert that the problem of 
pollution monitoring is not a sensor problem but a problem of mathemati- 
cal modeling and data processing . 

Part 1 of this report presents summaries of satellite sensor performance in 
three spectral bands, visible, infrared, and microwave. Part 2, the bulk 
of the report gives all the calculations, trade-offs and limitations of the 
three sensor systems. 
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PART 1: SUMMARY OF SATELLITE PASSIVE SENSORS 
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SATELLITE 


Visible 



Wavelength 
Spectral Band 
Sensor Number or Area 
Dwell Time 
Collecting Aperture 
Sensitivity 
Efficiency 
Spatial Resolution 
A1 titude 

Swath Width crosstrack 
downtrack 


510 nm 6B0 pm 

30 nm 
100 cm 2 
38s 

10 cm diameter 
Min. Det. Cont. = .04% 
1 % 

30 in x 30 m 

1 Nn 

1 Mm 
10 km 


SENSORS 


Infrared 


Microwave 


10. 6 iim 3.8wm 

2|im .8 win 

100 sensors 
1 50 ws 

30 cm diameter 
10 mk 50 ink 

50% 

100in x 1 00m 240m x 240m 

1 Min 

1 Mm 
NA 


1 cm (30 GHz) 
3 GHz 
N sensors 
150 N wsec 
10 m diameter 
.5K 


1 km x 1 km 

1 Mm 

1 Mm 
NA 




PART 2: SENSORS TRADE-OFF 
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In the following we present the results of our trade-off between sensitiv- 
ity and resolution for three satellite based sensors 

• visible imaging radiometer 

• infrared scanning radiometer 
%■ microwave scanning radiometer 
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1.0 


OPTICAL RADIOMETRY 


We have studied satellite sensors to investigate the feasibility of see- 
ing low-contrast patterns in color on the sea surface. The feasibility 
divides into two questions. First, what existing or planned sattelites 
provide a capability suitable for a demonstration. The short answer 
is that NASA LandSat D scheduled for launch In 1983 should provide 3Z 
contrast in a 30 x 30 meter pixel, slightly smeared in the crosstrack 
direction by limited electronic response time. 

The second question is. What capability can be achieved in a satellite 
package optimized for pollution detection? In this case NASA's present 
radiometric methods are inappropriate because they have traded off con- 
trast and resolution in order to achieve absolute radiometric calibration. 
If we forego absolute measurements, it is possible to optimize for recog- 
nition of low-contrast patterns that .-epresent spilled pollutant. This 
leads to a sensor that performs more like a low-light-level television 
system. 

The following sections address the two questions. The first describes 
the radiometer aboard LandSat 0, and the second describes possibilities 
for a system dedicated to the pollution detection problem. 

1 .1 LandSat D: Thematic Mapper 

This satellite carries a radiometer that scans a swath 100 nautical miles 
wide on the earth's surface. In the visible range of interest are two 
color bands, a blue-green one from 450 to 520 nm and yellow-orange from 
520 to 600 nm. Each band uses an array of 16 silicon photodiodes as de- 
tectors. A mechanical scanner sweeps the crosstrack direction while or- 
bital motion scans along track. 

The relevant properties of the LandSat D radiometer appear in Table 1 
along with the symbols used to denote these quantities in the aquations 


TABLE 1: PROPERTIES OF LandSat 0 


Swath width: S * 135 km 

Altitude: H * 1 Mm 

Velocity, equivalent surface: V * 6.6 km/sec 

Detectors in array: N * 16 

Footprint (square pixel) size: 5 * 30 meters 

Dwell time: t 3 lOusec 

2 

Spectral radiance of the sea: L x * 0.035 W/sr/m /nm 

Filter pass, 3and 1: b * SO nm 

Telescope efficiency: n * 0.2 

Aperture diameter: d 3 40 cm 

Detector responsivity: R 3 0.6/volt 

-21 

Thermal noise, spectral density: kT 3 4 x 10 watt/Hz 

(T 3 290 kelvln) 

- 1 9 

Electron's charge: e 3 1.6 x 10 coulomb 

Detector/amplifier bandwidth: B 3 50 kHz 

g 

Noise resistance: R 3 1.0 x 10 ohms 

(feedback resistor for transimpedance amplifier) 
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that follow. The area coverage rate Is the product of velocity V by 
swath width S: 


A - VS - 1220 km 2 / sec, 

From this we derive the dwell time r that one of the 16 sensors inspects 
each pixel (o by 6) on the sea surface: 

t * N6 2 /A * lOusec 

Note that N6 2 is the footprint of the N sensors on the earth surface. 

The scene (sea and air above it) has a certain radiance L (watt/m /sterad' 
ian) that is focused on the detector where power p produces current I: 

p * L (efficiency) (area) (solid ana'e) 

2 2 2 

. (solid angle) * (aperture area/altitude ) 3 -d /4H 
p 3 Ln6 2 (idV4H 2 ) 

£■ * J n ( f) Z - 22.6 x 10-’ 2 » 2 (1) 

The light detector is a silicon photodiode that converts electric current 
with typical responsivity R 3 0.6 amps/watt. To facilitate calculations 
that follow, we use one factor K to convert sea radiance to detector cur- 
rent I: 


K2 U = p L = R L a 13,6 x 10 " 12 m2sr/volt (2) 

We shall estimate the noise current 1 in the detector and use K to con- 

n 

vert it to a minimum detectable radiance (MDL) on the sea surface: 


MDL = i n /K 


(3) 


S 


Two types of noise limit performance of the system, shot and thermal 
(Johnson), given by the following: 

(4a) 
(4b) 
(5) 

Here R is the load resistance across the photodiode, or the feedback 
resistance in the case of a circuit with a transimpedance amplifier. 
Normally thermal noise would limit the performance of a photodiode, 
but in this case NASA has used a newly developed sensor with a trans- 
impedance amplifier. The equivalent noise resistance is remarkably 
g 

high, R a 10 ohms. 


1 2 • 2eIB 
i 2 * 4kTB/R 


Substituting Eqs. 4a and 4b in 5, and using Eqs. 2 and 3 gives 
MDt » J (C 1 C 2 L) B 


( 6 ) 


where 

C 1 * 4kT/K 2 R C 2 - 2e/K (7) 

and 


(^B = 4.35 x 10" 3 (W/m 2 /sr) 2 

( 8 ) 

C 2 B - 1.18 x 10" 3 W/m 2 /sr 

We have plotted Eq. 6 in Figure 1 in the form MDl/L, which is the mini- 
mum detectable contrast. One point on this curve is particularly sig- 

2 

nificant at a radiance of L $ = 2.3 W/m /sr. This value is typical of 
scene radiance looking at open ocean (no clouds). The minimum detectable 
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Minimum detectable contrast, percent 


3 



Fig. 1 . Minimum detectable contrast in 
thematic mapper, landSat D 
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contrast 3.1*, and Its reciprocal Is the signal-to-noise ratio in 

terms of photocurrent (I/i^), namely 32. It happens t h at NASA has 

specified this same signal-to-noise ratio at the same radiance in its 

contract with Hughes Santa 3arbara Research Center. The saturation 

2 

radiance is 10 W/m /$r at wnlch the specified signal-to-noise is 35, 
a bit higher than our estimate of 79 (reciprocal of 1.3* in Figure 1). 

We estimated the radiance of the sea by first averaging Its spectral 
radiance (see in Figures 12, 13 and 14 of Roswell Austin's paper in 
Optical Aspects of Oceanoorachv , edited by Jerlov and Nielsen and 

reproduced in the Appendix). The average spectral radiance L. * 

2 ' 

0.035 W/m /sr/nm, multiplied by the filter pass band, b 3 SO nm, then 

gives the result we used: 

2 

L $ a Lb a 2.S W/m /sr. (in the blue-green) 

The same gr?nhs indicate that only 142 (a seventh) of the light comes 
from the ■ surface; the rest is backscattered by air. Of the part 
from tne sea, about a third is surface reflection. Thus 102 of the 
radiance at the radiometer is light that has penetrated the sea. Pre- 
sumably this percentage is improved by avoiding sun glint. 

Finally, the frequency response of the thematic mapper rolls off at 
B * 50kHz * 1/2 t, and so the radiance in one pixel spills over a bit 
into the next one or two pixels in the crosstrack scan direction. 

1 .2 Dedicated Imaging Sensor 

This section discusses a satellite package dedicated to those applica- 
tions in which the main problem is to see low-contrast patterns. The 
best approach is to use a large array of sensors to increase the dwell 
time on each pixel (compared to LandSat D). An ordinary television 
camera has in effect an array of aoout 250,000 sensors occupying a total 

area of about one square inch. It would be reasonable to assume a total 

2 

photosensitive area as larce as 100 cm , which can be divided into vari- 


1 i 


ous numbers of pixels (TV lines) depending on the trade-off between con- 
trast and resolution. If common television camera tubes were used, 17 
of them would be required to achieve this area in each color band l.e. 

51 for a tri-color system. However, oversized tubes or some other ad- 
vanced solid-state image sensor would reduce the total number of devices. 
We shall leave the question of size open and denote the photosensitive 
area by a. 

A wide swath is very desirable to eaten as many breaks in the clouds as 
possible. We assume that the width equals the altitude, i.e. 

X * H * 1 Mm (9) 

as shown in Figure 2. In the direction along track, the field of view 
is smaller to avoid excessive photosensitive area, 

a = vy (10) 


(see Figure 2). If Y denotes the dimension of the field of view along 
track, then the dwell time is 

T = Y/V (Vi) 

It will be advantageous to make this time as long as possible, even 
10 seconds or more, to average out clutter from the sea surface. This 
feature requires some form of motion compensation, but the result is 
rewarding because it will very thoroughly remove statistical fluctua- 
tions in sun glint, and to some extent tho^e due to patches of foam. 
Cepending on sea state and wind, this can be a great advantage over 
previous NASA radiometers that essentially "freeze" the instantaneous 
clutter in each pixel'. 

Let us assume a fairly advanced system with fast optics, say a focal 
ratio of two: 
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f/d » 2. 


( 12 ) 


Any faster optics would not work well with an interference filter to 
limit the band to q.\ * b. Equation 12 determines the solid angle of 
rays into one point on the focal plane: 

Q * (tr/4) (d/f) 2 * >r/16 (13) 


Apart from an optical efficiency factor n Q . the image has the same radi- 
ance as the object, namely 

l s » L x b (14) 

where L is the spectral radiance and b the filter bandwidth. Thus the 
number of photons collected in the entire field of view is 


N * = (L s /hv)n o ai2 T 

and the number in each pixel is 

n . V 2 . 1 - n o L s ai ! 

% ' XY " T6 hv VH 


where Equations 9, 11, 13 and 15 have been used. 


(15) 


(16) 


The mean number of detection events in each pixel of the photosensitive 
surface is 


n , 3 n „ n . 
d q (j 


(17) 


where n q is quantum efficiency, and the standard deviation in this num- 
ber is 


An d 
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(18) 


These equations combine to give the minirum detectable contrast 


an. 


MDC 


► n n 
d < 


( 19 ) 


Substituting Eq. 14 in 16 in 19 gives 


MDC * 


,16 

\ * 



n_n_ 

O q 



1 

i*a 


( 20 ) 


Curiously, the collecting aperture d has dropped out of this expression 
because it is tied to f through the focal ratio, Eq. 12, and to other 
factors through the geometry of Figure 2. 


For plotting this expression, it is convenient to put 


L - £ A (21) 

A 

where A is the peak spectral radiance when the sky is clear *n'i the sun 
is high, 

2 

A 3 .035 W/m /sr/nm. 


and i represents all environmental losses of illumination due to sun 
angle, hydrometeors, or whatever. Finally we lump all losses of light 
together, whether in the environment or the equipment and put 

n * in 0 n d * (22) 

For example, we might have 

optical transparency - .5 (mostly in the filter) 

20% aperture observation - .8 
detector quantum efficiency - .15 


sun angle - .5 

haze and Rayleigh scattering - .33 


product: n * .01 


In this form, Eq. (2) becomes 


MDC * 


/16 (hv) VH 
' it nAb 


1 

3 > a 


( 23 ) 


Figure 3 shews MDC versus 5*a (5 in meters and *a in cm) with n as a 
parameter. 


For example, when 

2 

a 3 100 cm (10 by 10 detector surface), 
5 3 30 m, n * 15 

then MDC 3 3.3 x ID -4 . 


Quite clearly this is less than the natural contrast noise on the sea due 
to patches of foam, seaweed, windrows, and the like. Therefore, we desire 
a long dwell time on each resolution cell to average out these factors as 
much as possible. Equation 11 gives 

T 3 Y/V 3 (y/f) (H/V) (24) 

where we have used 


Y/H 3 y/f (25) 

which results from similar triangles in Figure 2. The factor y/f in Eq. 

24 cannot in practice be made arbitrarily large to attain a long dwell 
time because the motion compensation becomes nonlinear and complicated 
as the angles involved grow large. It is beyond the scope of this study 
to solve such details, but for the sake of illustration, suppose 
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Minimum Detectable Contrast 





y/f » 1/4, 


(26) 


which keeps the angles quite sma 1 1 (< 7° * arctan y/2f). Then Eq. 24 
gives T * 38 sec, a nice long period that averages out nearlv all foam 
statistics as well as sun glint. Moreover, using x * f (Figure 2) and 
xy * a (Eq. 10), we find that 

f - 2/a 

2 

If a 3 100 cm (previous example) then 

f 3 x 3 20 cm 
y * a/x * 5 cm 
d 3 f/2 3 10 cm, 

all very small and reasonable. 

For the case of an aircraft platform, we have not plotted the MOC as in 
Figure 3 because it is so small that it is meaningless, i.e. the product 
VH in Eqs. 20 and 23 is even smaller. However, the dwell time for aver- 
aging out clutter is not as long. Assuming 

H 3 3 km (9840 ft) 

V 3 160 m/s 

Equations 26 and 24 give T = 5 sec, long enough to average out sun glint, 
but not much else. 
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2.0 


INFRARED SCANNING RADIOMETRY 


For our purposes, "Infrared" means wavelengths from 3 to 20_m. Radiation 
from sea and clouds in this Pand consists almost entirely of thermal emis- 
sion, reflected and scattered sunlight being negligible. The so-called 
"photographic" infrared, wavelengths near l_m, should be classed with visi- 
ble light for our purposes, because this band is a component of sunlight, 
and so radiometry works only during the day. Within the band of interest, 
there are two particularly clear "windows", one at 10. Sum and the other 
at 3.8um, as shown in Figure 4. The 10. Sum band (9.6 to 11.6) is the bet- 
ter for three reasons. First, it provides 200 times as many photons, which 
overpower noise in the sensor and provide more sensitivity. Second, the 
long wavelength has minimum sensitivity to small particles, haze and such, 

4 

because Rayleigh scattering decreases as X . Finally, this window has the 
least atmospheric absorption (from gasses and vapors) of any optical band 
(including visible and ultraviolet). 

When the radiometer looks vertically down at the sea , it provides an almost 
pure measure of the temperature of the sea surface uncorrupted by other fac- 
tors. Reflected sky radiance comprises only 1 or 2 % of the signal at inci- 
dence angle of 0 to 30° , but this increases to 100% at grazing incidence; 
see Figure 5. The observed radiance originates very close to the surface, 
the mean distances being 

;(3.8) 3 60um 3 16x 

;(10.6) 3 Hum 3 lx (27) 

This means that the radiometer is very sensitive to any floating pollutant, 
such as an oil spill that inhibits evaporation at the surface. At the very 
least, a coating changes the temperature by 0.5K, which is the temperature 
chance through the so-called conduction layer , approximately the top 1.5 mm 
of the sea. This is a relatively large signal, but of course it is diluted 
with the signal from normal water wherever the film is thin and breaking up. 
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Fin. 4. Transmission over 







A sensitivity goal of one fiftieth this amount, l.e. 10 mK is probably ap- 
propriate. We show that an advanced satellite radiometer at 10.6um can see 
a (ICO m) resolution cell with this temperature sensitivity providin' t 
has 100 detectors In array. A radiometer at 3.8’.m would only resolve d..ut 
(240 m) . An airborne radiometer of the same design would be pvt- kill, l.e 
orders of magnitude excessive performance capability. 

We doubt if the Infrared radiometer has much sensitivity to pollutants that 
mix into bulk water, for then the surface effect is very small. There may 
be exceptions when the sun is shining on the water and the pollutant dis- 
colors it enough to raise the temperature by absorbing sunlight. Also, 
most chemicals release or absorb heat upon dilution, but this heat Is a 
transient occurrence observable only when the dump is quite fresh. 

2.1 Blackbody Formulas 

The sea Is very nearly a blackbody in the infrared with emissivity 

* 1 - * -98 to .99 

o o 

at normal incidence as discussed in connection with Figure 5. Thus it is 
appropriate to express radiance in terms of the blackbody function 3: 

L ( A ) * eB(x,T) power/area/steradian 

However, for reasons discussed in Section 2.3, it is more appropriate to ex 
press blackbody radiation in terms of photons instead of the tradiv.ona'i 
formulas for power. We denote this change by using the symbol Q instead of 
B: 


eB -*■ eQ, photons/sec/area/steradian 
We also denote spectral radiance with a subscript 

Q x * Q/b, ph/sec/area/sr/um (23) 
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where b Is the spectral bandwidth in micrometers. 


In these terms, the Planck radiation formula is 

* t! exp <• TT> 


(29) 


where 


3 . 2 . 

C 2 2c - 6.0E25 photons • -m /sec/m /sr 

(30) 

C g 2 hc/k * 143S8i.m • kelvin 

Quantities useful for evaluating Equations 29, 29, and 30 in the two bands 
of interest appear in Table 2 above the dashed line. Those quantities in- 
volving temperature are evaluated at 295 kelvin. The choice of bandwidth 
b (fifth in the list) is governed by the need to avoid absorption lines in 
the air spectrum. Figure 4. 

In all infrared radiometers, it is traditional to express the observed radi- 
ance in temperature units, i.e. the temperature of a blackbody that would 
give the same radiance as the observed value. This is quite natural because 
radiometers are calibrated by looking at a blackbody at a known temperature. 
Thus we need a formula to express a radiance increment aQ in terms f the 
equivalent temperature Increment: 

AT » aQ/Q' (31) 


where Q* 
sion: 


3Q/3T. Differentiating Equation 29 gives the required expres- 
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q 


TT ,n \ 



(32) 


or with the aid of Equation 31. 


TABLE 2 USEFUL RADIOMETRIC QUANTITIES 


evaluated a: two wavelengths and T » 295 Kelvin 


Quantity 

Units 

3.3un 

10.5'_m 

CJ\ • hv/k 

e 

kelvin 

3786 

1357 

CJ\T - hv/kT 
e 

(dimensionless) 

12.83 

4.601 

C q /\ 4 - 2c/A 4 

2 . 

photons/s/m /sr/um 

2.S8E24 

4.75E22 

Q x (Eq. 29 ) 

2 

ph/sec/m /sr/um 

7.67E18 

4.77E20 

b (suitable in air) 

urn 

0.6 

2.0 

Q(295,A) • Q x b 

2 

ph/sec/m 'sr 

4.6E18 

9.5E20 

Q/Q' * kT 2 /hv 

kel vin 

23.0 

64.1 

R a (nomal Incidence) 

AT 

m 

2.40 

0.83 

R fl Q/Q* 

millikelvin 

552 

532 


/ 

— note only 4U 
di fference 
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(33) 


(f )-f <T7» 

or 

f • f <r > ■ f <&> (34) 

e 

The ratio Q/Q * for Equation 33 Is listed below the dashed line In Table 2. 

2.2 Natural Background Fluctuations 

The last item in Taole 2 gives a particularly interesting special case, 
t.’te maximum effect that sky can have or. the apparent temperature of the sea. 
Suppose the sky is completely overcast, anu a radiometer beneath the over- 
cast looks directly down at the sea. (or else a satellite radiometer looks 
through a small hole in the overcast) Also suppose that the temperature of 
the overcast is not too different from that of the sea. (A small difference, 
say 5K, does not matter much since the sky contributes only 1 or 2* of the 
radiance.) Then the sea and sky in effect form an isothermal enclosure, 
and the radiometer necessarily sees the thermodynamic temperature of bo.h: 

^cloudy “ 

3ut when the sky is clear 

<>cle»r * £<) < T > " 0 ' R o> «< T > 

" ^cloudy " ^clear * 
and finally, using Equation 31, 

aT a $ 3 R o $*** 540 mK (35) 
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as listed in Table 2. If the sky is partly cloudy, then the temperature 
increment is 


dT(f) - f(R 0 Q/Q') (26) 

where f Is the fraction of cloud cover. It is not easy to determine f 
accurately. Since the radiometer sees a blurred image of the sky in the 
wind-roughened sea, one cannot say exactly how much of the *ky it is see- 
ing nor to what degree. If an estimate of f is in error 10°,, then the cor - 
responding error in temperature is 

AT - 0.1 R 0 Q/Q* * 50 mK, (37) 

which is about the maximum natural noise level for broken clouds , unless 
an auxiliary sky radiometer is used to make a correction. 

If radiometric temperature is measured at both 3.8 and 10.6um, then there 

are two Equations 36 whi-.h in principle be solved for two unknowns, water 

temperature and fractional cloud cover. But a curious coincidence occurs 

as shown in Figure 6. Tne x and y axes represent the two unknowns, and 

the locus of possible values is plotted for each radiometric measurement. 

The intersection that fixes f and T . has such an extremely small angle 

sea 

that the solution is worthless in the presence of expected noise. 

If the two radiometric temperatures, T(3.8) and T ( 1 0 . 6 ) , are substracted, 
then both of the principal effects, i.e. T_ and T. b .., are eliminated and 

' Sea SKy 

some other residual affect remains. We do not know what the principal con- 
tribution is, but it is not roughness, as our computer simulations have 
shown. It may turn out that the residual is very sensitive to pollution. 
Perhaps it is the temperature difference resulting from different skin 
depths in the two bands, Equation 27. The thermal gradient from evapora- 
tion and heat conduction at the sea surface is T * 3T/3z = .33 K/m, which 
gives 
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Fig. 6. Attempt to Find Both Cloud Cover and Water Temperature using 
Radiometric Temperature at Two Wavelengths 


27 



A g T(3.8) « 20 mK 


( 38 ) 


A g T(10.6) 3 4 mK ) 

If the sea is not too calm, say sea state 2 or greater, then the fluctua- 
tions due to clouds, aT * 50 mK, Equation 37 occur over larg» distances, 
comparable to the altitude of the cloud because the cloud image in the sea 
is blurred over about a radian. This may not interfere with the observa- 
tion of pollution plumes on a smaller scale, say a few hundred meters. 

Over these distances the natural background noise has a variance spectrum 
given very roughly by 

l(3E-3 cycles/m) 3 3E-3 K^tT ) 

(39) 

«(0.3 cycles/meter) 3 2.3E-3 K'm’" I 

(linear interpolation (ex rapolation) may be used for other frequencies). 
These values are very rough because they were derived from a 1 -dimensional 
cut through space time (a ship moving at a fixed speed) in a constant field 
of view, and so we could only assume that the fluctuations were frozen in 
the sea surface as the ship moved through. In a spacial bandwidth, Av x by 
AVy, the temperature variance is 

aT^ 3 $Av Av (40) 

a y 

2 

If we want the temperature variance in a (10 meter) resolution cell, it 
is reasonable to take 

* lv y * Ti T in ete ' rs * - 05 cycles/meter ('ll) 

Then Equation 40 gives 

AT 3 y 6E-3 lOrrT x .05 m’^ 3 4 mK (42) 

a much reduced natural background. 
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Aadiome te r Sensitivity 


? 3 

«■ • w 


A satellite radiometer must be quite sophisticated compared to an aircraft 
radiometer. The satellite sensor collects relatively few photons from each 
resolution cell on the sea surface because of its high altitude and short 
dwell time. We assume that the satellite radiometer employs low-noise am- 
plifiers of the type used in LandSat D and discussed in Section 1.1. In 
this case, the sensitivity is limited by shot noise (photon statistics) 
which overpowers other sources of noise, mainly thermal noise in the ampli- 
fier. (Considerably more photons are received in the infrared *han a simi- 
lar radiometer for visible light.) We assume that the infrared radiometer 
employs cold stops and a cold filter so that radiation from warm parts with 
in the radiometer cannot reach the detector. Thus the only shot noise is 
that from the field of view. 


The performance of infrared detectors and radiometers has traditionally been 

♦ 

expressed in power units using figures of merit such as NEP, D and D . How- 
ever, since the shot noise limit is both feasible and required, it is approp 
riate to bypass these quantities completely and discuss sensitivity in terms 
of photon statistics alone. It is for this reason that the preceding formu- 
las were expressed in terms of photon flux instead of power commencing with 
Equation 28. 


A straightforward range equation gives the number of photons collected from 
each resolution cell on the sea surface: 


n 

♦ 


n 


Q(\) 


/ d wel l \ / cel 1 \ / solid angle subtended 
\time / \area / l by radiometer aperture 


Q(t) (i 2 ) 



(«) 


where d is aperture diameter and H the radiometer's altitude. The number 
of detection events pet all is given by 
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n ■ nn 


(44) 

wnere n is an overall efficiency that includes optical transparency , aper- 
ture obscuration, and detector efficiency, the last being the efficiency 
in converting photons to charge carriers. 

Swell tima ? depends on N^, the number of detectors in array; the more there 
are, the more time each detector can spend cn one cell. Assuming an effici- 
ent scan pattern, i.e. each detector is always looking at a cell within the 
swath of interest, the dwell time is 

area examined at any instant 
area scan rate 


T 



(45) 


where S denotes swath width and V the platform velocity. Substituting Equa 
tion 45 into 43 into 44 gives 


n 


e 


, "N d Q 

4 SV 



(46) 


2 

Note the extreme sensitivity to resolution 5. This is because £ enters 
into both the emission area and the dwell time. 


Receiver sensitivity derives from the well known formula for the variance 
of photon (Poisson) statistics: 

(4n e ) 2 • n e (47) 

Converting to radiometric temperature uncertainty IT gives 


An e _ 1 . aQ . aT ( C e\ 

n e ‘ * n e “ Q “ T \'T/ 


(48) 
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wnere Equation 47 was used in the first step, 46 in the second, and Equation 
33 in the last step. Solving for aT gives 






(49) 


aT vT 4 S V 

T C t nTH? 

c u 


H 



(50) 


We have plotted this equation in Figure 31 in the form aT versus 5 with 
N, as a parameter. Even though aT is not very sensitive to we chose 
it as a variable because the number of detectors varies so greatly in di f 
ferent arrays, all the way from 1 to 200 in a recent NASA development. 
Reasonable values assumed for other quantities are the following: 


T 3 295 kelvin 
C e /\T, Q - See Table 2 
S * H * 1 Mm 
V 3 6.6 km/sec 
n 3 0.5 
d 3 30 cm 


Figure 7 shows that a 100-element array at 1Q.6um (solid line) provides 100- 
meter resolution and 10 mK sensitivity . At 3. Sum (dashed lines) the perform- 
ance is marginal, and either sensitivity, resolution, or both have to be 
compromised. On the left side of the graph the dashed lines are discontinued 
where the assumption of shot-noise-limited performance begins to fail. 
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Minimum Detectable Temperature Change, b.n 



3.0 


MICROWAVE RADIOMETRY 


We foresee a limited application of microwave radiometry for pollution sens- 
ing from satellites, although there may be some application from aircraft 
providing they fly beneath the clouds. The reasoning behind this conclusion 
follows. First, the reouirement for angular resolution forces the choice 
of a short wavelength on the order of \ » 1 cm (F * 30 GHz). Sut the radiome- 
ter is sensitive to clouds at this wavelength, and so its use on a satellite 
Is limited to areas with clear skys. Moreover, at these high frequencies, 
the radiometer is quite insensitive to sea temperature and conductivity (sal- 
inity), as discussed later. The main sensitivity that remains is surface 
roughness which is of great interest for detection and monitoring of oil 
si icks. 

However, the measurement of roughness puts radiometry into direct competi- 
tion with radar scatte. ometry , which measures the same thing. Radar has 
most of the advantages: 

1 . penetrates clouds 

2. high resolution available by range gating and synthetic 

aperture 

3. no problem with reflected cloud radiance 

4. lower frequencies (e.g. L-band) usable if synthetic 

aperture used for angular resolution 

5. cross polarization data 

The only advantages for radiometry are: 

6. no transmitter power 

7. no clutter statistics (random interference effects) 

In case of a satellite, Items 1 and 2 appear to be the driving factors. 

The cloud penetration (1) is very important because 70 to 30% of the ocean 
is covered by clouds. Also resolution (Item 2) is critical in the satellite 
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case as we shew In Section 3.3. Item 3, problems with reflected cloud radi- 
ance is severe because reflectivity is on the orcer of 50 f j, depending on 
frequency and look angle. This may tip the balance in favor of radar even 
in the aircraft case. The fact that radar requires transmitter power and 
has problems with clutter statistics (implied by Items 6 and 7) is just some- 
thing we shall have lO work around. These matters are described in the fol- 
lowing sections and summarised briefly in Table 3, a sort of 4-dimensional 
microwave matrix that considers: 

• platform: satellite, aircraft 

* mode: active (radar), passive (radiometer) 

e sky: clear, cloudy 

e sensitivity to water condition: roughness, salinity, 

temperature 

3 . 1 Radiometric Temperature of Seawater 

The Intensity of radiation recovered in a microwave radiometer is custom- 
arily expressed as radiometric temperature T. This would be the same as 
thermodynamic temperature if the radiometer were pointed at a perfect ab- 
sorber, but seawater is not, and so a radiometer looking at the sea sees 
a complicated mix of water and sky radiance reflected from various direc- 
tions. The results in the case of a clear sky is that the radiometric 
temperature depends very little on the water temperature except at low 
frequencies ; Stogryn [IEEE Transactions AP-15, p. 278 (Mar 67)] confirms 
it in a classic paper, for F =» 19.4 GHz (\ s 1.5 cm), windspeeds 0 and 
7 m/sec and nadir angles 0 and 50°, as shown in Figure 8. Other investi- 
gators have found similar results, e.g. a group at Texas A&M, Figures 9 
10, and 11. Figure 9 shows considerable dependence on sea temperature at 
9.3 GHz (X = 3.2 cm), but as discussed later, this frequency provides too 
little angular resolutions. c igures 9, 10, and 11 also show that the 
radiometer does not respond effectively to salinity, which might be import- 
ant in some pollution studies. 


Sensitivity 

water j conductivity I surface 
temperature) (salinity) | roughness 


Platform 


Satellite 


Ai rcraft 



Sky 

clear 

cloudy 

POOR compared to radar 
which offers far 
better resolution and 
the following options: 

• range gate 

• synthetic aperture 

• cross polarization 

NO Clouds 
require 
X £ 10cm 
but angular 
resolution 
requires 
X i 1 cm 

NO Sensitive only If 
X > 3cm, but angular 
resolution requires 
X 5 1cm. 

NO Sensitive only If 
X - 10cm, but angular 
resolution requires 
X ^ 1cm. 




OJ 

t/> 1 o 



clear 


YES but not 
superior to radar. 


cloudy (fly under?) 


POOR compared to radar 
since radiometer data 
complicated by reflected 
cloud radiance 


POOR since measure of conductivity 
must be disentangled from large 
effects of surface roughness and 
clouds 


Probably NO since Infrared radiometer 
does the job with far better resolution 
and less Interference from reflected 
sky and surface roughness. 


clear 


cloudy 


roughness 

YES 

conductivity 

YES If long X 

temperature 

NO 
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Ffg. 3. Temperature of horizontally 
polarized radiation as a 
function of water temperature. 
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Fig. 9. Theoretically pre- 
dicted curves. 
Brightness tempera- 
ture vs. salinity 
and temperature. 
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Fig. 10. Antenna brightness tempera 
ture of fresh water as a 
function of thermometric 
temperature. 



Fig. 11. Antenna brightness tempera- 
ture of sea water as a 
function of thermcmetric 
temperature. 
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A microwave radiometer looking at t.ne sea is sensitive :o tne ancle of the 
water surface, and hence to wind that roucrens the surface . This effect 
appears In curves by Stogryn reproduced here as Figures 12, 13, and 14. 

3.2 Clouds 

Figures 12 and 14 show that temperature sensitivity on the order of 
AT 3 1 Kelvin 

is needed to sense small chances in surface roughness. This in turn limits 
the amount of Interference that can be tolerated from clouds, aT c> This 
interference is proportional to the fractional absorption A: 

iT cloud ■ AT c 

Assuming T. 3 300 Kelvin gives A 3 .003. Comparing to Figure 15 shows that 
a light cloud causes appreciable interference at frequencies 10 GHz 
or higher. But frequencies this high (see Section 3.3) are required to 
attain reasonable angular resolution, and so the microwave radiometer, un- 
like radar, simply will not be very useful in seeing though clouds. 

* i 
» 

3.3 Instrumental Limitations 

The sensitivity of a microwave radiometer to temperature changes is given 
by 

2 T n 

AT * - (51) 

• 

where T n is the receiver noise temoerature, 3 its bandwidth, and t the 
dwell time. The factor of 2 applies strictly to a Dicke radiometer, and 
varies slightly with other schemes used to stabilize the gain. The dwell 
time needed for Equation 1 depends on platform velocity V and geometrical 
factors: 


I 

I 
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ATTENUATION NEPERS 



Fig. 14. Temperature of 
norizontally polarized 
radiation as a function 
or wind speed. 


**Uiy 


Or 



Fig. 15. Zenith microwave opacity components for a standard atmosphere, 
3 gm/cm^ columnar water vapor density and two clouds. 
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resolution are a 
area searen rate 


t - 5 2 /SV (52) 

where S is the swath width and 5 the dimension of a resolution cell (pre- 
sumed square). At nadir, the resolution is approximately: 

5 * Xh/d , (53) 

where h is altitude and d aperture diameter. This equation expresses a 
severe bind: X should be short for good resolution, but long to avoid 

excessive sensitivity to thin clouds or fog. The best one can do is use 
large apertures as listed in Table 4 along with assumed X, h and resulting 
5. The value 6 * 1 km in the satellite column is adequate for mapping a 
large oil spill, but too large for many smaller dumps. Other reasonable 
assumptions are listed in this table along with the results from Equations 
51, 52, and 53. These results show that the sensitivity is quite adequate. 
The real crunch is in resolution and sensitivity to clouds. 




TABLE 4 RADIOMETER 

ASSUMPTIONS AND 

RESULTS 

Assumptions 


Satel 1 ite 

Aircraft 

Altitude H 


1 Mm 

3 km 

Swath S » H 


1 Mm 

3 km 

Velocity V 


7 km/ sec 

160 m/s 

Wavelength \ 


1 cm 

3 cm 

Frequency F 


30 GHz 

10 GHz 

B * F/10 


3 GHz 

1 GHz 

Aperture diameter 

j 

10m * 1000 A 

1 meter * 33\ 

Noise temperature 

T n 

1 50K 

1 50K 

Results 

Resolution 6 


1 km 

100 m 

Dwelling time x 


140us 

20 msec 

Sensitivity AT 


0.46K 

.067K 
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APPENDIX 

Curves of Spectral Radiance 
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Figures by RW Austin, "Remote Sensing of S 
Aspec ts of Oceanography, edited 
Academic Press, NY 1974 (continued next 
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